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Abstract 

We present geometric methods for uniformly discretizing the continuous A''-qubit Hilbert 
space Hm- When considered as the vertices of a geometrical figure, the resuhing states form 
the equivalent of a Platonic solid. The discretization technique inherently describes a class 
of 7r/2 rotations that connect neighboring states in the set, i.e. that leave the geometrical 
figures invariant. These rotations are shown to generate the Clifford group, a general group of 
discrete transformations on A*' qubits. Discretizing Hm allows us to define its digital quantum 
information content, and we show that this information content grows as N'^. While we believe 
the discrete sets are interesting because they allow extra-classical behavior — such as quantum 
entanglement and quantum parallelism — to be explored while circumventing the continuity of 
Hilbert space, we also show how they may be a useful tool for problems in traditional quantum 
computation. We describe in detail the discrete sets for one and two qubits. 

1 Introduction 

The discrete nature of the configuration space for N classical bits is the key property allowing 
robustness of digital computation. The Hilbert space for N qubits, on the other hand, is a 
continuous complex manifold. This continuity appears essential to the exponential speed-up of 
some quantum computing algorithms, such as Shor's factoring algorithm[l], over their classical 
counterparts^. But it also poses a challenging problem for the experimentalist: errors in quantum 
gates are themselves continuous, so even minute errors can accumulate throughout the execution 
of an algorithm and lead to its failure. 

Yet, quantum error correction and fault-tolerant computation schemes have been developed to 
meet this challenge[2, 3, 4]. That reliable quantum computation is possible using both a noisy 
quantum register and noisy gates is a result of surpassing importance. However, such schemes 
still place stringent fidelity requirements on the basic quantum gates and the quantum register: 
estimates for the threshold error probability above which they fail are typically 10~^-10~^[5]. 

Can universal control of a scalable quantum register with this level of fidelity be realized? If 
so, are there concepts we can borrow from digital computation that might facilitate the develop- 
ment of this technology? If not, are there intermediate computational paradigms that might relax 

^As evinced by the Gottesman-Knill theorem, for example. See [11]. 
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these requirements, but still exploit "extra-classical" phenomena such as quantum parallelism and 
quantum entanglement? We note that these arc still possible in a discrete Hilbert space. 

Quantum gates are typically implemented by applying time-dependent fields to a qubit system. 
They correspond to rotations of a unit vector in Hn, with the angle of rotation usually determined 
by the duration and amplitude of the pulse which generates the field. In principle, such rotations 
are simple to implement, given an appropriate time-dependent Hamiltonian. But in practice, noise 
in both the qubit system and applied fields inevitably leads to errors. Sophisticated techniques that 
build up a desired gate from a sequence of rotations about successively orthogonal axes have been 
developed to mitigate the effects of noise. In the field of NMR, especially, techniques for performing 
high-fidelity rotations are now very mature [6]. Yet such techniques for protecting against noise are 
not directly gcneralizable to arbitrary angles and axes of rotation. As a result, most experimental 
protocols for quantum manipulations rely as heavily as possible on a small set of rotations, usually by 
an angle of 7r/2 or tt, specifically optimized for the given qubit system. In the landmark experiment 
by Vandersypcn, et al., in which an NMR-bascd quantum processor was used to factor the number 
fifteen via Shor's algorithm, the protocol contained a single rotation by an angle less than 7r/2 — a 
conditional 7r/4 rotation [7]. 

Nonetheless, \miversal control of a quantum register requires in theory only a finite number of 
discrete gates, provided the gates form a universal set. Then an arbitrary "software- level" quan- 
tum gate can be constructed to a precision e by concatenating 0(log'^(l/e)) discrete gates from 
the imiversal set (c w 2) [5]. However, given that each discrete gate itself would likely comprise 
a sophisticated series of rotations, the prospect of concatenating C'(log^(l/e)) such gates to cre- 
ate each software-level operation — and doing so before the register decoheres — makes the fidelity 
requirements of fault-tolerant computing schemes all the more exacting. 

Much of this difficulty in achieving high-fidelity control of a quantum register can be alleviated 
by limiting ourselves to nonuniversal sets of quantum gates which generate only finite transfor- 
mation groups. A finite transformation group implies a finite number of possible states, so this is 
equivalent to imposing a discretization on the underlying Hilbert space: the quantum register be- 
comes "digital" . By suitably choosing the transformation group, the allowed states can be selected 
to have certain well-defined properties, such as known expectation values with respect to a set of 
measurement operators. 

As an illustration of this idea, consider the task of testing the experimental protocol for gen- 
erating one-qubit rotations, which can be represented on the Bloch sphere. Suppose we wish to 
optimize the fidelity of a 7r/2 rotation about the y-axis in a given qubit system. Starting from the 
state |0), we perform a counter-clockwise 7r/2 rotation about the y-axis, yielding the target state 
= (|0) -f \l))/^/2, then we perform a measurement in the {|0) , |1)} basis. By repeating this 
many times, we obtain the expectation value {az) of the target state |-t-a;). Assuming imperfections 
in the state preparation and readout have been accounted for, this expectation value would ap- 
proach zero if our 7r/2 gate were perfect, since (+a;|(T2| + x) = 0, while deviations from zero would 
imply an imperfect 7r/2 gate. Specifically, (az) = S would imply that, on average, the gate has per- 
formed a rotation by an angle of 2 arccos \/i6 + l)/2. With knowledge of other one-qubit discrete 
states and their expectation values ((7;) , we could also test rotations about the x- and z-axes. By 
direct generalization, this simple protocol can be used to test rotations on any number of qubits, 
provided we have an appropriate discrete set of target states. 

Also, in the nascent field of quantum feedback control, techniques have been developed to 
dynamically correct quantum processes. By performing continuous weak measurements on the 
quantum system, it is possible to control and correct quantum state evolution through feedback[8, 9]. 
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Incorporating these techniques into quantum computing experiments could also be facilitated if the 
number of processes involved were reduced to include only a small class of rotations connecting 
states with well-defined properties. 

Having explained why we wish to consider a discrete subset of the full continuous Hilbert space, 
we would now like to draw a geometrical analogy. Since we do not want to privilege any region 
of Hilbert space over any other, the set must comprise a uniform sampling of H^. The structure 
of the finite sets we have in mind is exemplified in real space by the Platonic solids — geometrical 
figures such as the tetrahedron, cube and octahedron characterized by the geometric equivalence 
of their vertices — which represent discrete subsets uniformly spanning a sphere in M?. In short, we 
are seeking to generalize the Platonic solids to Hilbert space by selecting from Hj^ a finite set of 
states corresponding to the vertices of a 2^-dimensional complex uniform polj^tope. We call such 
subsets uniform Hilbertian polytopes, and denote them by ^ jv . 

In discretizing an A/'-qubit register, what extra-classical phenomena must be sacrificed? If we 
select the transformation group to include only it rotations — the quantum generalizations of the 
NOT gate — we generate only a discrete set of 2^ states, and fall back on a purely classical register, 
with no possibility for quantum entanglement. But, as we will show, the next level towards finer 
rotations, the transformation group based on 7r/2 rotations, is sufficient for rich extra-classical 
behavior: the number of discrete states in the set then grows as 

2(N -H3JV)/2[io], the majority of 

which are entangled for N > 2. Also, the super-extensive growth of the discrete set relative to the 
classical nuinb(;r of states 2^ implies that a great deal of the quantum parallelism possible in the full 
Hilbert space remains possible in the discrete set. Though such a digital quantum register would not 
allow algorithms which are exponentially faster than their classical counterparts (Gottesman-Knill 
theorem), a possible reduction of an algorithmic scaling speed from 0{N'^) to 0{N) could still be 
useful. 

At the same time, within the framework of traditional quantum computation, a discrete set and 
its associated transformation group can provide a useful "roadmap" for navigation in the entire 
Hilbert space. 

The notion of discrete sets of AT-qubit states is not novel. Indeed, discrete sets have already 
been considered in quantum error-correcting codes. There, a special set of 2* orthogonal states, 
to be used as codewords for the basis states of k encoded qubits, are selected from a higher- 
dimensional continuous space H^- Gottesman's stabilizer formalism provides a general framework 
for describing and producing quantum error-correcting codes, and allows an analysis of a broad 
class of quantum networks in the Heisenberg picture[ll, 12]. Powerful though it is, the stabilizer 
formalism approaches the problem of discretization algebraically; it does not address the geometric 
relationship between the discrete quantum states, nor the relationships among the various gates 
that connect these states. 

The purpose of this paper is thus to provide such a geometric approach to the uniform dis- 
cretization of i^AT, and to suggest the use of such discrete sets, either as an arena for exploring 
extra-classical behavior, or as a heuristic tool for the analysis of certain quantum information pro- 
cessing problems. We refer to these notions collectively as digital quantum information^ . 

For simplicity, we focus here on the one- and two-qubit Hilbert spaces. However, most of our 
results are directly generalizable to higher-dimensional spaces. When possible, we use a language 
that makes this generalization straightforward, if tedious. In section 2, we treat the discretization 
problem using stabilizer theory and derive a class of generalized 7r/2 rotations belonging to the 

^Later, we will use the phrase rigorously, in reference to the information content of the discrete Hilbert space. 
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Clifford group that connect states in the discrete set. Later in section 3, we present an alternate, 
purely geometric approach to discretization based on shelling the high dimensional lattices. 



2 Discretization based on stabilizer theory 
2.1 Stabilizers and the generalized Pauli matrices 

We begin this section with some essential results from stabilizer theory. First, define the iV-qubit 
Pauli group Gn as the set of all Ai'-fold tensor products of 2 x 2 Pauli matrices, with four possible 
overall phases to satisfy the closure requirement: 



where^ 



{(Tw,ax,ay,az}'^ (g){±l,±i}. 



1 

1 

1 

1 
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Clearly, each element of Qn acts on the A/'-qubit Hilbert space. Qn has order 4^+^, and is generated 
by a minimal set of 2N elements, i.e. two non-identity cr's acting on each qubit. We refer to 
individual elements of C/at as generalized Pauli matrices, and denote them Sc«/3...c = (ya®(^i3®' • ■®(yc,- 
The generalized Pauli matrices share many of the properties of the 2x2 Pauli matrices. For example, 
they all either commute or anti-commute, and 

Sj = Sj (Hermitian), 

= id (Square root of unity), 
TrStSfc = 2^5jfe (Orthogonal). 

A stabilizer is an Abelian subgroup of the Pauli group. In the present work, we are predominantly 
concerned with the commutation properties of the generalized Pauli matrices, so we neglect the 
phases {±1, ±z} required for closure of Qn under multiplication. That is, we deal with the set Sn 
of 4^ generalized Pauli matrices rathe;! than the group Q^- To distinguish the Abelian subsets of 
iSjv from the Abelian subgroups of ^jv, we refer to the former as psewdostabilizers, a name which 
also highlights the close relationship between this work and stabilizer theory. 

The largest possible subsets of Sm whose elements all mutually commute have 2^ elements. 
These maximal pseudostabilizers will form the foundation of our first discretization procedure. 



2.2 The uniform Hilbertian polytope 9y 



N 



We are now in a position to discuss a formal definition for the uniform Hilbertian polytope for N 
qubits. First, we establish the desired properties the discrete sets must have. We seek to construct 
S^N, such that: 

''The present numbering scheme has been chosen to coincide with the binary vector space representation of 
stabilizer codes, as in [12]. 
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1. It contains all the states |6o&i-"^iv-i) corresponding to the classical bit configurations. 

2. Each state of ^^at is geometrically equivalent to all the others (uniformity). 

3. The distance between two normalized states and defined as 



satisfies^ 



djk > 7r/2 for all i,k. 



4. It is the largest set of states which satisfies the above requirements. 
Denote by s% the maximal pseudostabilizers in .Sjv- Then these desired properties are obtained 
if we adopt the following construction for the vertices oi Sjn'- 

Definition 1 An N-qubit state vector is an element of^jpf if and only if it is a common eigenvector 
of each element of a maximal pseudo stabilizer . That is, if is a generalized Pauli matrix on 
N qubits belonging to s^, is an N-qubit state vector, and Xj is an eigenvalue ofT,j belonging 
to the vector , 

1*,) e^jN^ 1*,) = A, 1^-,) for all e s^. 

As a consequence of this definition, and from the theory of stabilizers, we find: 

a) Each s^, which has 2^ — 1 elements different from the identity, generates 2^ different 
discrete states all separated by djk = tt. Each state corresponds to a unique pattern of Xj = ±1. 

b) Each shares exactly half, or 2^~^, of its elements with its nearest neighbors; 2^~^ 
with its second-nearest neighbors, etc. Any discrete state in ijjv therefore has N "levels" of non- 
orthogonal neighboring states. 

c) For each and each of its nearest neighbors one can associate by a general algorithm 
a transformation from the common eigenvectors of to those of s^. That is, any two states of 
Sjn are linked by a finite sequence of similarity transformations. 

d) The similarity transformations are formed from generalized orthogonal tt/2 rotations of 
the form: 

X^i = i= (Sfe + iS() , where Sfc, e s%. 
v2 

The superscript a denotes a subset to which both its S's belong and the subscripts specify the 
S's. The inverse operations are: 

This definition implies that for any X, 



X^X = \A (Unitary), 
X* = -id (7r/2 Rotations), 

which is consistent with the property that a spin- 1/2 acquires an overall phase of e™ = — 1 when 
rotated by 27r. 

^Discretizations with a finer minimum distanee may be useful and would be interesting to explore (for two qubits, 
see section 3). For one qubit this could correspond, for instance, to the icosahedral geometry. 
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e) The X's generate the Chfford group Cn, defined as the normahzer of the Pauli group[12], 
which has the property of leaving ^ n invariant (proof to follow) . 

f) The set Sn of generalized Pauli matrices on N qubits contains 

N-l 

s = II (2^-'= + 1) 
fe=o 

maximal pseudostabilizers s^. Each has 2^ elements, and contributes 2^ simultaneous eigenvec- 
tors. The uniform Hilbertian polytope on N qubits Sj n therefore contains 

iV-l 

= 2^ JJ (2^-'= + 1) 
fe=o 

vertices, or states [10]. The following table gives the first values of Vn, along with the number of 
classical bit configurations for comparison. 



A' 


1 


2 


■■'> 


4 


~) 


(i 


i 


Vn 


6 


60 


1080 


36, 720 


2,423,520 


315,057,600 


81,284,860,800 


Cn 


2 


4 


8 


16 


32 


64 


128 



The digital quantum information in N qubits can be defined as the information content of Sjn, 
i.e. as log2 Vjv. It is easy to show that Vn grows as 2^^ +3iv)/2^ ^^lis information content is 
super-extensive in N. While it is insufficient for algorithms which would be exponentially faster 
than classical ones, it is nonetheless a remarkable property for a discrete space. 

We now turn to an explicit construction of the uniform Hilbertian polytope for the one- and 
two-qubit cases. 



2.3 The one-qubit case, S^i 

We show here that Sji is isomorphic to an octahedron. For one qubit, the set Sn is simply the Pauli 
matrices: <Si = {ayj,aj:,ay,az}- The last three a's anti-commute with one another, while they all 
commute with the identity cr^. So the three sets of mutually commuting matrices are trivial to 
construct: si = {aw,az} , = {aw,aoc} and sf = {a^, Oy} . 

When the elements of s\ are diagonalized, we obtain the computational basis: 

l+^> = |0), 
|-^) = |1)- 

s\ generates the pair 

|0) + |1) 



\+x) = 
\-x) = 



V2 ■ 

|o>-|i> 



x/2 

And s\ generates 
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\+y} = 
\-y) = 



|o) + ^|i) 

\o)-^W 
V2 ■ 



There are three orthogonal 7r/2 rotations, which form the seed of S^i: 

^01 = ^9 (^o+icTi), 



V2 

^03 = TsK + ^fTs)- 



The diagonalization of the seed elements leads directly to the six eigenstates, as summarized in the 
table below. The states are listed here as unnormalized row vectors for clarity, and are separated 
into columns according to their eigenvalues. 



Set 


1 + i 


1 - i 


1 


(1,0) 


(0,1) 


2 


(1,1) 


(1,-1) 


3 


(1,^) 


(l,-») 



Each of the 7r/2 rotations has an inverse: 



(^0 fe) ' = -«^fe " 71 ^"^^ ~ ^"""^ ' 



It is easy to verify that 



-{(7o+l(7k) = l(Jk 



and that the X's are mapped into one another by similarity transformation: 

1 1 



^0 j^Oi 



((To + icTj) ((To + iUi) ((To - iUj 

(2cro + 2i eijkCTk) 



2V2 
1 

2V2 



■^a k if ^ijk 



1 



= (^Ofe) ^ if ^ijk = -1 



This implies that each X transforms a member of iji into its neighbor. 
Proof. // 

and if 
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then 

= A, 1^,0-)) 

therefore is an eigenvector of X'^ . ■ 

The X's with their inverse generate a twenty- four element group isomorphic to the octahedral 
group of pure rotations which leaves the octahedron invariant. 

2.4 The two-qubit case, 5^2 

The set of generalized Pauli matrices for two qubits ^2 comprises 4^ = 16, 2^ x 2^ matrices given by 
^Aai = cx® o'h, \yfi = w, X, y, z, as presented below. We write the index on the a's in binary, then 
concatenate the two strings to form the new index for the E's. For example, cr^ cr^ = era cti = 
CTii (8) ctqi — Siioi = S13. 
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So 



^10 



Si2 



Syx — Si 



1 








1 
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1 














1 
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-1 
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i 














—i 





— i 
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i 









^wz — Si 



Sujy — S3 



Sxz — S5 



Szj/ — S7 



Sx2 — Sg 



Sxy — Sii 



Sj/2 = Si3 



Sj/y — S15 
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—i 
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The products of these matrices can easily be found from 
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The maximal pseudostabilizers in S2 are presented below. ^ 





Letter notation 


Number notation 


1 

1 


/V V V V 1. 


iy„ V, y, y^x 


9 


ry V V V T. 

X'^wwi '^wx^ ^zw: ^zx) 


/v„ v„ y. y^.x 


Q 


Jv V V V T. 


5y„ y,-, y, y^x 


I 


SV' V' V V V 


f V , V , V V 1 
1^0: ^1: ^1)/ 




SY y y y \ 
X^wwt ^wx: ^xw: ^xx) 


■[Zjo, Zj2, Zjg, ZjioI 


U 


iv V y y \ 

X^wwi ^wy; ^xw7 ^xy J 




7 


\^ww^ ^wz^ ^ywt ^yz} 


{So, El, E12, E13} 


8 


^^WW! ^wxt ^yw: ^yx} 


{Eo, E2, E12, E14} 


9 


\^ww: ^wy: "^yw: Eyy} 


{Eo, E3, E12, E15} 


10* 


\^WW7 ^zz: ^xx: ^yy} 


{Eo, E5, Eio, E15} 


11* 


\^ww: ^zz: ^xy: ^yx} 


{Eo, E5. Ell, E14} 


12* 


\'^WW7 ^zxj "^xz: Eyy} 


{Eo, Eg, Eg, E15} 


13* 


{E^^, ^ZXj Ea;y, Syjg} 


{Eo, Eg, Ell, E13} 


14* 


{E^^, Yj^y^ Y^xz: '^yx} 


{Eo, E7, Eg, E14} 


15* 


\^WWi ^zyi ^XXi Ey^} 


{Eo, E7, Eio, E13} 



Each of these fifteen sets, or pseudostabilizers, will yield four simultaneous eigenvectors, contributing 
four states to ^2- We therefore recover the result that SJ2 has sixty states. 

These subsets can be classified as corresponding to entangled or product states by examining 
their generators. Note that each pscudostabilizcr is generated by any two of its non-identity ele- 
ments. The presence of the one-qubit identity aw when the generators are decomposed into tensor 
products of one-qubit Pauli matrices implies that the states corresponding to that subset are prod- 
uct states. Conversely, the absence of the identity in this decomposition indicates that the states 
corresponding to that Abelian subset are fully entangled states. The subsets whose generators do 
not contain the one-qubit identity are denoted here and below by an asterisk(*). 

We can obtain all the states of Sj2 directly by forming a mixed linear combination of the first 
two non-identity elements from within each set. For instance ^3^12 = ^ (E^y + iEy^,), when 
diagonalized, gives four orthogonal eigenvectors with four different eigenvalues. (Note that other 
rotations from si, such as Xi2,i5 = (E12 -|- iSis) and Xs^is = ^ (E3 -|- iEis) will produce the 
same four eigenvectors, but with permuted eigenvalues.) We thus construct in this manner fifteen 
generalized X's, each having a different principal axis, which form the seed of ^32- Diagonalization 
of the seed X's exhaustively gives the eigenvectors constituting S}2- We list these eigenvectors 
below, separated into columns corresponding to the eigenvalues (±1 ± z). For clarity, we list them 

^Constructing the subsets of mutually commuting Pauli matrices can be done through a series of logical steps. 
The key is to note that [Sjfc, ^im] = requires either [iTj, c;] = [cfc, <Tm] = or {<7j,ai} = {ufc, Cm} = 
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Figure 1: Graph of the set ^2 of two-qubit gcncrahzcd Pauh matrices Ej (circles bearing the 
subscript of the matrix in letter notation) and by the pseudostabilizers (triangles formed by three 
connected circles). The three E's in a triangle share four common eigenvectors which form an 
orthonormal basis spanning the two-qubit Hilbcrt space. The fifteen triangles thus give fifteen sets 
of four basis vectors. Shaded triangles correspond to entangled states while non-shaded triangles 
correspond to product states. Neighbouring triangles have one (non-identity) S in common, and 
each (non-identity) S is shared by three triangles. The line segments joining the vertices of a 
triangle correspond to pairs {j,k} of commuting matrices; each segment therefore specifies a 7r/2 
rotation Xj^k = (Sj + zSfc)/\/2 that transforms the eigenvectors of an adjacent triangle into its 
neighbour. The figure thus constitutes a "roadmap" for navigating the discrete set Sj2- Repeated 
circles indicate the closure of the graph. 
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as unnormalized row vectors. Again, entangled states are denoted by an asterisk on the set label. 





— 1 — t 


1 J- i 


1 ■} 


1 -1- i 


1 
1 


{0, U, U, I) 


[U, 1, U, U) 


l^U, U, 1, 


1^1, U, U, Uj 


9 


[U, u, — i, ij 


( — L, i, U, Uj 


(^U, U, i, ij 


U, J^, u, uj 


Q 
O 


u . u , /. , 1 J 


[ I.. 1 . u , U I 


(n n i ^ \ 


— /- , 1 . u , U J 


A 


((1 1 (1 1 "1 

[ u , — 1 , U . 1 ] 




i" 1 (11 n"! 




O 


(,-'-5 J-i J-i J-j 




('—I _1 1 \\ 
I J-I J-I J-, 


('1111') 


u 






—1 — 7 

I', J-, *, J-J 




7 


(0,i,0,l) 


(0,-i,0,l) 


(t, 0,1,0) 


(-Z, 0,1,0) 


8 


(-«,z,-l,l) 


(«,-*, -1,1) 


(j,i,l,l) 


(-Z, -2,1,1) 


9 


(-l,z,z,l) 


(1, -1,1,1) 


(1, «,-«,!) 


(l,z,z,-l) 


10* 


(0,-1,1,0) 


(-1,0,0,1) 


(1,0,0,1) 


(0,1,1,0) 


11* 


(z, 0,0,1) 


(0, -i,l,0) 


(0,2,1,0) 


(-2,0,0,1) 


12* 


(1,1,-1,1) 


(-1,1,1,1) 


(1,-1,1,1) 


(1,1,1,-1) 


13* 


(*,-*, 1,1) 


(«,i,-l,l) 


1,-1) 


(-2,2,1,1) 


14* 


(i,l,-z,l) 


(-i,l,z, 1) 


(z,-l,z,l) 


(2,1,2,-1) 


15* 


(-l,-i,i, 1) 




(l,-i,-i,l) 


(1,2,2,1) 



Note that the first stabilizer of the product sector corresponds to the computational basis, while 
the first stabilizer of the entangled sector corresponds to the Bell basis. 

This method not only finds the states of in an exhaustive way. It also provides a road map 
for navigating the discrete set. To illustrate this, consider three pscudostabilizcrs which wc call s°, 
and s°. They have two generalized Pauli matrices in common, one of them being the trivial Eq. 
We chose one of the two and call it S„. Consider three S ^ S^, 

e s", Efc G s^ and G s'' 

such that 

{Sj, Efe} = {Efe, S;} = {El, Tij} = 0. 

Since all pairs of generalized Pauli matrices that do not commute must anti-commute, they also 
satisfy 

Then it is easy to show by a direct calculation that 



^™i = ^(S™+iS,-), 



have one of the properties: 
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Y'f^ Y'tt 

mk m j km ml 

or 

va Y^b _ / Y^c \-l _ 
mk mj km \ ml) Im 

depending on whether the two anti-commuting cr's in the decomposition of SjSfe appear in cychc 
order or anti-cycHc order, respectively. This means, following the proof given for the one-qubit case, 
that all the eigenvectors of .s" arc transformed into the eigenvectors of s'^ by the transformation . 
All together, there are 120 different generalized two-qubit 7r/2 rotations generated by the scheme 

Xi j = — ^ (S, + iSj) , 

where Sj and commute. 

Among these X's, there is a subset that plays an important practical role. These are 7r/2 
rotations of the form 



Xo j ^ —r= (So + iT,j) . 



1 

They correspond to the unitary time-evolution operator 



U {t) = e'^-""" 

with T ~ 7r/4, and are thus directly implemented by a Hamiltonian proportional to Ej. These 
rotations constitute the practical means of navigating Sj2- They can be seen as the "primitives" of 
the Clifford group, as we show below. 

But first, it is important to note that the generalized Pauli matrices in the above arguments arc 
not limited to the two-qubit case, but can in fact be over any number of qubits. These results are 
therefore directly generalizable to larger Hilbert spaces Hn and larger discrete sets f)N- We duly 
conclude that our generalized 7r/2 rotations on A'' qubits, constructed from the pseudostabilizers 
s^, leave Sjn invariant. 



2.5 The generalized tt/2 rotations generate the A^^-qubit Clifford group 

So far we have successfully discretized the continuous Hilbert space Hj^, and in doing so we have 

described a class of generalized 7r/2 rotations that leave the invariant. From the point of view of 
operators acting in Hn, this discretization means we have reduced the continuous transformation 
group SU{2^) to a finite group. Here we show that this finite group is the A'-qubit Clifford group 
Cjv. 

The Clifford group is defined as the normalizer of the Pauli group. That is, a unitary operator 
X is contained in Cjv if and only if 

First, let us show that our X's are elements of Cat. That is, 

Xj k = A={^j + i^k) e Cn if P,-, Sfc] = 0. 
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We have 



where sji = ±1 and eui = il.Thus, 



If e/cjeji = 1, 



If Ekieji 



2 

-SjiT,i {T,j + iekiSji'S'k) {^j — i'^k) ■ 



= -Eji'^i (2So + iSfcSj - iSjSfe) 



= —lejiTii'SjYlk € ^jv- 

So the generaUzed 7r/2 rotations on A/' qubits are elements of the Chfford group. 
Now note that the Chfford group is generated by the Hadamard, 



'1 1 
1 -1 



phase, 
and CNOT, 



CNOT 



+ i 
l-i 
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1 
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gates [5]. The Hadamard gate may be composed from the one-qubit 7r/2 rotations X02 and Xqi: 

H = Xq2Xqi{X2o) ^ 

= —iXo2XoiXo2 
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The phase gate may be trivially constructed from a single one-qubit 7r/2 rotation: 



-'^01 

1 + i 






1-i 



while the CNOT is simply the product of three X's : 

CNOT = {Xo2)-^Xo e{Xo 4)-^ 
= —X20X0 6X40 
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So our generalized 7r/2 rotations allow a direct construction of a gate set that generates the Clifford 

group. The finite transformation group leaving 9)n invariant, generated by the generalized 7r/2 
rotations on N qubits, is thus the A^-qubit Clifford group Cjv- 



2.6 Comments 

One of the motives we presented for this work was the difficulty we anticipate in achieving the 

reliability requisite for fault-tolerant quantum computation. Clearly, limiting the register to a finite 
number of possible states must alleviate this difficulty, but by how much? 

It can be derived from the properties of the Pauli group that each pseudostabilizer has 
N levels of non-orthogonal neighbors. Since the eigenstates of neighboring pseudostabilizers are 
connected by a single 7r/2 rotation, any state on can be reached from any other in at most A^-|- 1 
such rotations. This is to be compared with the result that an arbitrary state in the full Hilbert 
space can be reached to within an error e by concatenating C'(log'^(l/e)) rotations from a universal 
set, with c « 2. In addition, note that the CNOT and H gates are not directly implemented by a 
physical Hamiltonian, but must be built up from 7r/2 rotations which are naturally realized with 
accessible field variations, so there is a second simplification from working with the 7r/2 rotations 
rather than standard universal gate sets such as {H, S, T, CNOT}. 

From an experimentalist's point of view, therefore, the X's form a very natural language for 
building quantum gates. A rotation of the form X^j = (Eo + i'^j) I \/2 is directly implemented by a 
term in the Hamiltonian proportional to S^-. And as shown above, this class of rotations generates 
Cat. The X' s are thus the basic instructions for a sort of "machine language" for quantum processors. 
The following section shows a simple example of their calculus. 
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2.7 Sample application of digital quantum information: implementing 
CNOT 

The Hamiltonian describing a given physical system determines which of the generaUzed 7r/2 rota- 
tions will be directly realizable in that system. Implementing CNOT according to the decomposi- 
tion in section 2.5 requires a physical system with a Hamiltonian proportional to Y,zx in order to 
realize the entanghng operation Xq q = (T^ww +i^zx)/V^- Though this type of inter-qubit interac- 
tion is possible^, most qubit systems rely on a less exotic interaction, such as one proportional to 
Hzz or T^xx- How can we implement the CNOT gate in one of these more standard registers? 
Specifically, suppose the system is described by a two-qubit Hamiltonian of the form 

H = a{t)T.i + 6(i)S2 + c{t)T,4 + d(t)S8 + e(i)S5, 

where the tuning parameters a, b, c, d, e allow the relative strengths of the terms to be adjusted dur- 
ing an experiment. Our task is to replace the rotation Xq q in the sequence {Xo2)~^Xo e{Xo 4)"^ 
with a rotation or sequence of rotations generated by the above Hamiltonian. Following the dis- 
cussion in sections 2.3 and 2.4 on the relationships between the generalized tt/2 rotations, it is 
straightforward to calculate that Xqq = {Xq 3)~^Xo 5X0 3 while Xq 3 = {Xq i)~^Xo 2X0 1. To- 
gether these give an alternate decomposition that employs only directly realizable X's : 

CNOT={Xo2r^Xoe{Xo4)-^ 

= (-^02)~"^(Xo 3)~^Xo 5X0 3{Xo 4)""^ 

= (-'^O 2)~^Xo 1X0 2{Xo l)~^Xo 5{Xo i)~^Xo 2X0 l{Xo 4:)~^- 

Since the decompositions differ in their implementation but not in their meaning, such se- 
quences arc "synonyms." It is important to note that the CNOT is not a special case: a synonym 
suitable for a particular implementation could likewise be calculated for any gate in C2 . 

2.8 Conclusion to section 2 

In this section we have presented a geometric method for producing from the continuous Hilbert 
space Hn a discrete, uniform sampling S)n. Because all the states in the discrete set are geometri- 
cally equivalent, represents a generalized Platonic solid in Hpf. This method is closely related 
to the stabilizer formalism of quantum error-correcting codes. Inherent in our construction is a 
description of how different elements of jv are related by transformations generated by physical 
Hamiltonians expressed in the basis of generalized Pauli matrices. This has been demonstrated in 
detail for ^1 and ^2, and is obtainable by direct analogy for higher-dimensional spaces. These ideas 
provide a useful tool for analyzing problems in traditional quantum computation, as the example 
above illustrates. And though computation over the discrete set i3Ar is clearly less powerful than 
computation in the full Hn, it is potentially more powerful than classical computation. 

3 An alternate approach to discretization: shelling the high- 
dimensional dense lattices 

In this last section, we present an alternate approach to discretization that addresses the Hilbert 
space directly, without reference to operators in or the relevant transformation groups. 
^As a charge-flux coupling between superconducting qubits, for example. 
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Our strategy here is the following (we note n = 2^): the normalization condition, together 

with the writing of complex numbers as pairs of real numbers, identifies the Hilbcrt space Hn to 
the high-dimensional sphere S'^n-i embedded in R^". In order to discretize these hyperspheres, we 
use the successive shells of dense lattices in M?". At the same time, we must take into account the 
global phase freedom, and show how a discretization of the projective Hilbcrt space is induced. This 
means that several points on S'^n-i ^jjj represent the same physical state, as explained below. In 
light of this, it is important to distinguish between "qubit states" — the quantum states associated 
to the points on S'^"~^ and "physical states" -the states in the projective Hilbert space, which 
has the geometry of a complex projective space CP"~^. 

As in section 2, we again focus on the one- and two-qubit cases. The discretization of Hi'm first 
presented in terms of the 24 vertices of a self-dual polytope on S^, denoted {3,4,3}, which is the 
first shell of the densest packing in M^, denoted A4. These 24 vertices corresponds to 24 one-qubit 
states, and, modulo a global phase, to 6 physical states. For the two-qubit case, we use the Gosset 
polytope, the first shell of the densest packing in R^, denoted Eg- We find that this polytope has 
240 vertices corresponding to 240 two-qubit states, which leads to 60 physical states: 36 separable 
states and 24 maximally entangled states, just as in section 2. 

The A^-qubit Hilbert space is high-dimensional, and its multipartite nature (it is the tensor 
product of single-qubit Hilbert space) induces a subtle structure related to the state's various levels 

of entanglement[13, 14], which is not fully understood for A > 3. The Hopf map that is used here 
in the two-qubit case is entanglement sensitive, which translates here in grouping sets of equally 
entangled qubit states. 

3.1 The one-qubit case 

The generic one-qubit state reads 

1*) =fo|0)+ti|l), 

with \to\^ + 1^1 1^ = l.The normalization condition identifies the set of normalized states to a sphere 
embedded in M^. The projective case — the set of states modulo a global phase — leads to the 
Bloch sphere description, which can be seen as the base of the Hopf fibration,[15, 16]. An 
interesting discrete model on is provided by the self-dual {3, 4, 3} polytope [17]. It is related to 
the "Hurwitz" quaternion group. We now give two possible (dual) coordinates for its vertices, in 
each case as a real quadruplet and a complex pair. The correspondence between real quadruplets 
and complex pairs amounts simply to taking the first two (last two) real numbers as the real and 
imaginary part of the first (second) complex number. The first (second) complex number in the 
pair corresponds to to (ii). 

A first set, denoted Ti, is the union of the eight permutations of type (±1, 0, 0, 0) and the sixteen 
permutations of type ^(±1, ±1, ±1, ±1). Note that, modulo a global phase factor, these twenty-four 
points really represent six different physical states, which appear on the Bloch sphere as opposite 
points on the three orthogonal axes x, y, z. Indeed, the four points. 



Real quadruplets 


Complex pairs 


(1,0,0,0) 


(1,0) 


(-1,0,0,0) 


(-1,0) 


(0,1,0,0) 


(z,0) 


(0,-1,0,0) 


(-z,0) 
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represent the states |^i,a;) = |0), with u = 0, 7r/2,7r, 37r/2 ,which map to the same point on 
the Bloch sphere (the north pole), and they are therefore associated to the physical state I'^i) . 
Equivalently, the four points 



Real quadruplets 


Complex pairs 


(0,0,1,0) 


(0,1) 


(0,0,-1,0) 


(0,-1) 


(0,0,0,1) 


(0,i) 


(0,0,0,-1) 


(0,-z) 



represent the four states \'^2,(^) = e*"|l) with uj = 0, 7r/2, tt, 37r/2. The other sixteen vertices 
represent four other physical states, in the following way: 

1*3) ^ (|0) - |1)) , 1*4) ^ ^^^9^ (|0) - |1)) , 

1*5) ^ (|0) + i |1)) , 1*6) ^ (|0) - I |1)) , 

with ti) = 0, 7r/2, tt, 37r/2. 

For the later purpose of a discrete two-qubit construction, it is useful to describe a second 
version of the polytope {3, 4, 3}, for which the twenty-four vertices form a set T2 given by twenty- 
four permutations of the type {±1, ±1, 0, 0}/a/2. This polytope is obtained from the former one 
through a screw motion on S^oi angle 7r/4. This set leads to twenty-four states 

|$;,w) =e|*;,w), l = 1..6, w = 0,7r/2,7r,37r/2, and e = e'''/^ 

and to the six one-qubit physical states |$;) identical to |*/) . Indeed, the six states |*j) sit at the 
vertices of a regular octahedron. Since the states I*;, w) only differ from |*(,a;) by a global phase, 
they map onto the same six points on the Bloch sphere. 

3.2 The two-qubit case 

We now consider the two-qubit case, for which a generic state reads: 

I*) = too |00) + toi |01) + tio |10) + tn 111) , and \toof + |ioi|' + \tiof + \tiif = 1. 

The normalization condition identifies the set of normalized states to the sphere 5^ embedded in 

As for the one-qubit case, we consider the first shell of points in the densest lattice in M*, denoted 
Eg. This lattice belongs to the family of laminated lattices A^, and is therefore sometimes denoted 
Ag. These laminated lattices form a series which starts with the triangular lattice in 2d (the densest 
lattice in 2d). A3 is obtained as a particular sequence of A2 lattices packed in a third dimension, 
which gives the face centered cubic lattice, one of the two densest lattices in 3d. Appropriately 
packing A3 lattices along a fourth dimension leads to A4, whose first shell is precisely the {3,4,3} 
polytope we used above. Upon iteration, this construction eventually leads to the Ag = Es lattice 
suitable for the two-qubit case. We shall focus here on the set of 240 sites belonging to the Eg first 
shell that forms the Gosset polytope and, as for the one-qubit case, enumerate the physical states 
they represent. 
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3.2.1 Discrete Hopf fibration for the Gosset polytope on S"^ 

The 240 vertices of the Gosset polytope belong to the sphere 5^. These 240 vertices may be 
separated into ten equivalent subsets, each belonging to non-intcrsccting spheres. This is nothing 
but a discrete version of the 5'' Hopf fibration, with fibers and base S'^flS, 16, 18, 19]. 

It is simpler to use here quaternionic coordinates instead of complex or real ones. The above 
set Ti, scaled such that the corresponding points belong to a sphere of radius now reads: 

where i, j and k are the standard unit quaternions. The set T2 stays on a unit sphere and reads: 
T, = {i=(±l ± i), i=(±l ± j), -^(±1 ± k), ^(±i ± j), ^(±i ± k), ^(±j ± k)}. 

The 240 vertices of the Gosset polytope belong to the ten sets: 

5i = (T2,0), 52 = (0,T2), 53 = (Ti,Ti), 54 = (Ti,-Ti), 55 = (Ti,iTi), 
56 = (Ti,-iTi), 57 = (Ti,jTi), 58 = (Ti,-jTi), 59 = (Ti,kTi), 5io = (Ti, -kTi). 

Each of the ten sets gives a copy of a {3, 4, 3} polytope on a fiber S^. The points can be Hopf 
mapped, as described elsewhere [13], onto the base space S^. The location of the mapped point 
is intimately related to the entanglement of the corresponding two-qubit state. Indeed, the Hopf 
map is simply described as a first map which sends the pair (gi, 52) onto the quaternion Q — qiq^^ 
(which is sent to infinity if (72 =0), followed by an inverse stereographic map which sends Q to 5*^. 
A main result is that the Hopf map is sensitive to entanglement: for separable states, Q is simply a 
complex number, not a generic quaternion; conversely, for maximally entangled states (MES), the 
purely complex part of Q vanishes. This translates onto the base S*"* in the following way. Embed 
S'^ into R^, with coordinates {x;, ^ = • • • 4}; then separable states are such that Xs = Xi = 0, and 
the S'^ sphere spanned by {xq, xi, 0:2} form the standard Bloch sphere of the first qubit. Maximally 
entangled_states map onto the unit circle in the plane (ss, 2:4). Note that a well known entanglement 
measure, the concurrence [20], is simply given by the radius in the plane {x^^Xi): c = v^cf+xf, an 
expression which will be used later. 

In the present case, it is then easy to verify that the sets to Sq correspond to separable 
states, while sets 6*7 to correspond to maximally entangled states. The correspondence between 
vertices and states is made by transforming back the quaternion pairs into complex quadruplets 
whose terms are {too,toi,tio,tn)- More precisely, the {qi,q2) pair reads {too + toij, ho + tiij) ■ 
Note that the quaternion unit j acts on the right of the complex numbers, while it acts on the left 
in the definition of ^7^8. Since quaternion multiplication is non-commutative, this distinction is 
important in going back and forth between the lattice points and the states. 

3.2.2 The separable states 

Consider the set ^i, corresponding to the twenty-four states such that tio = tu = 0, and which 
reads 

[0)1 O [$i,a;)2, w = 0,7r/2,7r,37r/2, and I = 1..6. 
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As a whole, the six sets Si - ■ ■ Sq encompass 6 x 24 = 144 vertices, forming altogether 36 physical 
states, with four values of the global phase for each qubit state. Note that the precise value of 
the phases are important here in order that our discretization procedure uniformly cover the full 
Hilbert space. Using the above defined eigenstates of the one-qubit Pauli matrices, these states 
read: 



|±a;) |±a;) e*W4+™V2) 


\±x) O \±y) e*W4+mV2) 


|±a;) (g \±z) e*™'^/^ 


\±y) ® |±a;) e*('^/4+m^/2) 


\±y) (g) \±y) e*(V4+mx/2) 


\±y) ® \±z) e"™'^/^ 


\±z) ® ±4 e*""^/^ 


\±z) ® \±y) e™'^/^ 


\±z) ® \±z) eH'r/4+m7r/:i) 



where m = 0, 1, 2, 3 triggers the global phase. Each of the nine entries stands for the four possible 
sign combinations, leading to the announced thirty-six physical states. A simple view of these 
separable states consists in relating them to the "product" of two octahedra, each one belonging to 
the Bloch sphere of the individual qubits. 



3.2.3 The maximally entangled states 

The remaining four sets (altogether 4 x 24 = 96 sites) lead to a slightly more subtle structure. 
Wc find a total of twenty-four different physical MES, with four phase-distinct two-qubit states for 
each. But in the present case, the phase-distinct states actually belong to two different sets, either 
{Sr,Ss) or {Sq,Sw) 

As an example, we consider the set Sr and enumerate the states corresponding to, say, the 
quaternion pair (^-^y Tsj) ^ ^■'^^'"■'^ translates into the complex quadruplet (:^)0,0, ■^) and 
therefore to the MES: 

\MESuO) = (|00) + 111)) = i= (1+^, +z) + \-z, -z)) . 

The shortened notation \+z, +z) stands for |+-2)j(S) 1+0)2 • There is only one other element in S^, the 
pair (^—-^, "Tfj)' corresponding to \MESi,Tr) = e*'^ \MESi,0) . The other two elements belong 
to the set Ss and read and (^^, , associated respectively to the states \MESi,Tr/2) 

and |M£^5i, — 7r/2). Before giving the full list of states, it is interesting to focus on the next three 
states in St, generated by i, j and k. One gets, respectively, 

\MES2, 0) = -J= (|00) - 111)) = i\+z, +z) - \-z, -z)) , 

\MES„ 0) = i= (|01) + |10)) = i= {\+z, -z) + \-z, z)) , 

\MESi,Q) = (|01) - |10)) = {\+z, -z) - 1-0, z)) . 

The set {\MESi,Q) ,\MES2,Q) ,\MES3,Q) ,\MES4„0)} forms the entangled "magic" basis de- 
scribed elsewhere[21, 22]. Each physical state \MESi) corresponds to the four two-qubit states 

\MESi,^) , with a; = 0,7T for states in Sj and uj = tt/2,—tt/2 for states in S^i. As a whole, the 
forty-eight vertices of the sets ^7 and Sg provide twelve physical MES, each one representing a set 
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of four phase-related states. They are listed below, in the w = case, and with e = e"^/*: 
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- (loo) - 


\^o)^ + 

|iu;j i- 
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|ii) 
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'•2^) — y)) ) 
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— Z, —X)) , 


\MESq,0) = 


|(|oo)- 


|oi)) + 


e 
2 


|10) + 


|ii) 


e 

^71* 


\+z, -x) - i 


-2;,a;)) , 


\MESio,0) = 


|(|00) + 


|10))- 


e 
2 


101) - 


|ii) 


e 

^71' 


\+z,+y) + \- 


-z,-y)), 


\MESii,0) = 


f(|00) + 


|io»- 


e 
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|oi)- 


|ii) 
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^71* 


l+^,-y) + |- 


z,y)) , 


\MESi2,0) = 


|(|oo)- 


|oi)) + 
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110) + 


|ii) 


e 
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The latter twelve physical states may be read off from the two sets and 5io: 
MES,3, 0) = |00) + -1= 111) = -L (1+^, +z) + i\-z, -z)) , 



MESi5,0) 
MESie,0) 
ME 8x7,0) 



;jj|01)-^|10) 



1 

1 
1 

71 



(1+^,-2;) -\\-z,+z)) , 
(i \+z,-z) - \-z,+z)) , 



2 

M£;5i8,0) = |( 

Mi;5i9,o) = |( 

M£;52o,0) = |( 
MES2i,0) = |( 



M£;6'22,0) 
Mi;S'23,0) 
M^524,0) 



00} + |01)) - - (|10) - 111)) = {\+z, +x) - \-z, -X)) , 
00) - i |10)) + I (|01) + i 111)) = {\+z, -y) -i\-z, +y)) , 
00) - i |10)) + I (|01) + i 111)) = {\+z, +y) - i 1-2, -y)) . 
00) + |01)) - i (|10) - 111)) = 1= {\+z, +x) + \-z, -X)) , 
00) - |01)) + I (|10) + 111)) = {\+z, -x) + \-z,x)) , 
00) + i |10)) - I (|01) - i 111)) = {\+z, +y) + i \-z, -y)) , 
00) +i |10))-l(|01)-i |ll)) = -^(|+z,-y) + i \-z,y)), 
00) - |01)) + i (|10) + 111)) = {\+z, -X) - \-z,x)) . 



Succinctly, the ninety-six two-qubit entangled states can be written: 

\MESi,u;}=e''^\MESi,0}, with a; = 0,7r/2,7r,37r/2 andZ = l---24. 
Modding out the global phase, we can write the twenty-four physical states in the form: 



^[\+z,+z) + e^'\-z,-z)) 


^^[\+z,-z) + e^''\-z,z)) 


^{\+z,+x)+e^'^\~z,-x)) 


^{\+z,-x)+e^''\-z,+x)) 


^{\+z,+y) + e^'^\-z,-y)) 


^{\+z,-y) + e^'^\-z,+y)) 



with = 0,7r/2,7r,37r/2. 

Note that these twenty-four entangled states, together with the above thirty-six separable states, 
are in one-to-one correspondence, up to a global phase, with the sixty discrete states on Sj2 presented 
in section 2. 
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3.3 Comments 



3.3.1 Finer discretizations of H2: higher Eg shells 

Thus far, this alternate technique has provided compelling confirmation of the results from section 2. 

However, the present purely geometric approach does not describe the discrete set's transformation 
group. This point could be in principle addressed, albeit in a much less transparent way than in 
the first part of this paper, by using the Eg lattice point group (which has a high order) and finding 
the subgroup of rotations which survive the Hopf map. But this lattice approach does have the 
benefit of allowing discrete sets with a finer minimum distance (i.e. < djk < n/2) to be explored 
in a straightforward manner: while the first shell of Eg provides the discrete set with a minimum 
distance of 7r/2, a finer discretization might be achieved by considering the higher order shells in 
Eg,. This construction would provide a uniform set of two-qubit states, some of which would have 
intermediate entanglement. A note of caution is in order here, since we are only interested in 
describing normalized quantum states. Lattice points which are aligned, as viewed from the origin, 
contribute to the same two-qubit state. One should therefore only focus on the "visible points," 
which form the lattice's Mobius transform[23] . 

We do not give here a detailed description of these finer discretizations of H2. However, we 
note that the number Mj of sites on the J*** shell around an Eg vertex is simply given by [24] 

Mj = 240 

d\J 

where d denotes integers which divide J. The table below displays these numbers for the first four 
shells. Again, the physical states are obtained from these two-qubit states by modding out a global 
phase. 



J 


1 


2 


3 


4 


M 


240 


2160 


6720 


17520 



The shell by shell analysis, and its relation to the Hopf map, was done elsewhere[16, 18]. It 

allows us to get points on the second shell corresponding to states having concurrence 0,1/2,1 /V^, 1. 
The third shell contributes states of concurrence 0, 1/3, 2/3, V5/3, \/8/3 and 1. 

3.3.2 Discrete one-qubit mixed states 

A second advantage of our lattice approach over the pseudostabilizer approach is that it allows a 
discussion of discrete sets of mixed states. It is well known that the full set of one-qubit mixed 
states can be obtained by tracing out one qubit of generic two-qubit pure states. Mixed states are 
represented by points inside the so-called "Bloch ball," bounded by the pure-state Bloch sphere. In 
the context of generalizing the Bloch sphere for two-qubit pure states using the S**" Hopf fibration[13], 
it was shown that the Bloch ball corresponds precisely to the restriction to the triplet (xq, xi^x^) on 
the base space 5*. This describes mixed states obtained upon tracing out the second qubit. With 
this in mind, we are tempted to propose, in parallel to the two-qubit pure state discretization, an 
Sg-related discretization of the Bloch ball. 

From the Eg first shell, one gets the six permutations (±1,0,0), corresponding to pure states 
on the Bloch sphere (the traced separable two-qubit states) forming an octahedron. But one also 
gets the point (0, 0, 0), the Bloch ball centre, corresponding to traced maximally entangled states. 
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Then, from the Es second shell, we find in addition the eight permutations (±i, ±i, ±i) — a cube 
of radius •\/3/2 — and the twelve permutations (±^,±|,0) — a cuboctahedron of radius l/-\/2. The 
Bloch ball can be further discretized by using traced states originating from higher Es shells. 

3.3.3 Conclusion to section 3 

We have presented a second technique for producing a uniform discrete set of states from the 
continuous Hilbcrt space, again focusing on the one- and two-qubit cases. The pseudostabilizer 
and dense lattice strategics were done independently, and their agreement provides a compelling 
confirmation of our results, while also calling for a better understanding of their possible deeper 
relationship. While the dense lattice strategy is less easily focused on the transformation groups 
leaving the discrete sets invariant, it does have the advantage of allowing a discussion of discrete sets 
with a smaller minimum distance between states, and also of discrete sets of mixed states. Looking 
ahead, a discretization of higher-dimensional Hilbert spaces may be achieved in the dense lattice 
approach by using high-dimensional lattices in 2^ dimensions [25]. Though it does not describe the 
entanglement properties as nicely as in the two-qubit case, the S^^ Hopf map, corresponding to the 
three-qubit case, has been found to be entanglement sensitive[26, 27]. This case should be related to 
the dense lattice Aig in M}^, and is presently under study. It is interesting to note that the number 
of lattice sites closest to the origin — the lattice "kissing number" — for this case is 4320, which is 
precisely four times the expected number of vertices on the uniform Hilbertian polytope .$33. We are 
therefore likely to face a similar situation as in the one- and two-qubit cases, where there were four 
phase-related qubit states associated with each physical state. However, the four-to-onc relation 
between the Aig first shell sites and the vertices of Sj^ remains to be checked. Generalization to 
more than three qubits cannot use the Hopf fibrations, limited to S^^. A particularly interesting 
family to be checked further is the one described long ago by John Leech[28], which coincides with 
those studied here for iV = 1, 2 and 3, and whose kissing number is, for any N, precisely four times 
that given in the first part of this paper for the number of states in the generic Hilbertian polytopes. 
As for the three-qubit case, this precise four-to-one relation remains to be checked. 

M.D. and C.R. are indebted to Ike Chuang for his enthusiastic teaching of quantum information 
theory and to Daniel Gottesman for communicating and explaining his results. In addition, R.M. 
would like to thank Perola Milm.an and Karol Zyczkowski for interesting discussions and comments 
on the two-qubit Hilbert space geometry, and Philippe Biane who, while reading the first version 
of this manuscript, drew his attention to John Leech's 1964 paper. Daniel Esteve and Denis Vion 
are also gratefully acknowledged by all of us for helpful interactions. Finally, John Preskill's web- 
accessible notes and exercises on quantum information have been very useful. This work has been 
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